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Abstract 



Motivated by questions in mass-action kinetics, we introduce the notion of vertexical 
family of differential inclusions. Defined on open hypercubes, these families are charac- 
terized by particular good behavior under projection maps. The motivating examples 
are certain families of reaction networks — including reversible, weakly reversible, endo- 
tactic, and strongly endotactic reaction networks — that give rise to vertexical families 
of mass-action differential inclusions. We prove that vertexical families are amenable to 
structural induction. Consequently, a trajectory of a vertexical family approaches the 
boundary if and only if either the trajectory approaches a vertex of the hypercubc, or a 
trajectory in a lower-dimensional member of the family approaches the boundary. With 
this technology, we make progress on the global attractor conjecture, a central open 
problem concerning mass-action kinetics systems. Additionally, we phrase mass-action 
kinetics as a functor on reaction networks with variable rates. 
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1 Introduction 



The global attractor conjecture has been a central open problem in reaction network theory 
since its formulation by Horn in 1974 [14] . It asserts that any complex-balanced mass-action 
kinetics system of ordinary differential equations with positive initial conditions possesses a 
globally attracting stationary point in each stoichiometric compatibility class (see Conjec- 
ture [6J] for a more precise statement). It is well-known that this conjecture is implied by 
Feinberg's persistence conjecture a version of which asserts the following: for weakly 
reversible networks taken with mass-action kinetics, no species asymptotically becomes ex- 
tinct or unbounded. An a priori special case of the persistence conjecture asserts that at 
least one species survives asymptotically. In this paper, we show that this a priori special 
case in fact implies the entire persistence conjecture (Corollary I6.2[) . As a consequence, if 
the persistence conjecture is false, then in every minimal counterexample each species be- 
comes either extinct or unbounded. It follows that the persistence conjecture in dimension 
n implies the global attractor conjecture in dimension n + 1. 

Pantea [18], building on earlier work by Craciun, Nazarov, and Pantea [7], used the 
persistence conjecture in two dimensions to prove the global attractor conjecture in three 
dimensions. Pantea's work relied in part on projecting trajectories to lower-dimensional 
faces. Here, we generalize this projection argument in two ways. First, our main re- 
sult, Theorem 13.131 applies in arbitrary dimensions. Second, our results hold not only for 
mass-action kinetics networks, but for certain families of differential inclusions that we call 
vertexical (Definition I3.12H . The notion of vertexical family makes precise the essential 
structure required of a family of dynamical systems on hypercubes of varying dimensions to 
permit a structural induction argument of this sort. Theorem 13.131 shows that a trajectory 
in a vertexical family approaches the boundary if and only if either it approaches a vertex 
of the hypercube, or a lower-dimensional trajectory in the family approaches the boundary. 

Vertexical families of differential inclusions arise naturally in reaction network theory by 
way of mass-action kinetics on networks that are reversible, weakly reversible, endotactic, 
strongly endotactic (Definition I4.9l4j) . and so on. We prove that these networks, and more 
generally, projective classes of networks (Definition 15. ip . give rise to vertexical families of 
mass-action differential inclusions (Theorem 15. 191 and Corollary l5.20p . In the course of prov- 
ing this result, we are led to view mass-action kinetics as a functor (Theorem 15. 17j) . Doing 
so clarifies the concept of vertexical family and suggests that other questions concerning 
mass-action kinetics systems may be amenable to structural induction. Finally, Section [6] 
discusses the implications of our results for persistence of mass-action kinetics systems. 
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cussions, and Duke University where many of the conversations occurred. 

2 Dynamical properties of differential inclusions 

In this section, we recall certain dynamical properties of differential inclusions defined on 
manifolds. For background on manifolds, see [17] . All manifolds considered here have finite 
dimension. 

Definition 2.1. Let M be a smooth manifold. Let 7Tm '■ TM — > M be the tangent bundle 
of M. A differential inclusion on M is a subset X C TM. 

Example 2.2. The simplest differential inclusions on M are vector fields on M. The subset 
X C TM for a given vector field is the image of the corresponding section M TM. 

Definition 2.3 (Trajectory). Fix a differential inclusion X on a smooth manifold M. Let 
I Q K>o be a nonempty interval (in particular, connected) containing its left endpoint. 
A differentiable curve / : I — > M is a trajectory of X if the tangent vectors to the curve 
lie in X. An unbounded interval is a ray. A trajectory defined on a ray eventually has a 
property P if there exists T > such that property P holds for the function for all t >T. 

A smooth manifold with corners M is a space locally modeled on the closed nonnegative 
orthant [17^ p. 363]. The relative interior M \ dM of a manifold with corners is a smooth 
manifold [HI p. 386, Ex. 14-19]. 

Definition 2.4. Fix a smooth manifold with corners M with relative interior M = M\dM, 
and let V C dM be a subset of the boundary. A differential inclusion X C TM is 

1. persistent relative to V if the closure in M of every trajectory of X is disjoint from V. 

2. repelled by V if for every open set 0\ C M with 7 C Oi, there exists a smaller 
open set O2 f= Oi with V C O2 such that for every trajectory / : / — >• M of X, if 
/(inf J) ^ 0\ then /(/) n O2 = 0; in other words, if the trajectory begins outside 
of 0\, then the trajectory never enters Oi- 

3. permanent relative to V if it is persistent and there is a compact subset Q C M\U such 
that for every ray I, every trajectory of X defined on / is eventually contained in £1. 

More generally, a set X of differential inclusions on M is persistent relative to V, repelled 
by V, or permanent relative to V if every member X € X has the corresponding property. 

Remark 2.5. Any differential inclusion repelled by V is also persistent relative to V. 
The converse is false in general because different trajectories starting outside 0\ could get 
arbitrarily close to V. However, if M is compact, then it can be shown that this does 
not happen without violating persistence. Therefore, repelled and persistent are the same 
property for compact manifolds with corners. 
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Remark 2.6. A differential inclusion that is permanent relative to V need not be repelled 
by V when M is not compact. (When M is compact, see Remark 12 .51 ) The reason appeared 
already in Remark 12.51 different trajectories starting outside 0\ could get arbitrarily close 
to V. The crux is that the eventuality condition defining permanence does not apply to 
trajectories defined on finite intervals. On the other hand, if every trajectory defined on a 
finite interval can be extended infinitely — that is, can be extended to a trajectory defined on 
a ray — then in fact V does repel any differential inclusion that is permanent relative to V. 

Conversely, a differential inclusion repelled by V need not be permanent relative to V , 
even if M is compact, because O2 might necessarily be smaller when 0\ is smaller. That is, 
trajectories that start closer to V could eventually remain closer to V, as in Example 12.71 

Example 2.7. Fix a differential inclusion X on a planar disk M whose trajectories form 
concentric circles about its center. X is repelled by the boundary circle V = dM. Indeed, 
if 0± C M is an open set containing V, then the compact set M\0\ achieves a maximum 
radius r from the center, so we can take O2 to be the set of all points in M of radius > r. 

Lemma 2.8. If the differential inclusion X in Definition \2.4\ is persistent, then for every 
trajectory f : I M of X, there exist disjoint open sets Of and Oy in M containing the 
closures in M of f(I) and V respectively. 

Proof. A manifold with corners is metrizable, and hence it is a normal Hausdorff space. □ 

Remark 2.9. Even with Lemma 12.8^ our definition of persistence is weaker than the stan- 
dard definition |12] that requires each coordinate of a trajectory to remain bounded away 
from for all time. However, we follow other works in the chemical reaction network liter- 
ature [3j [5] in using Definition 12.41 1 . Additionally, Definition 12.41 1 agrees with the sense in 
which Feinberg first stated the persistence conjecture (see Section [6]). 

Mathematically, the definition we adopt has the advantage of being purely topological, 
so it behaves well under homeomorphism. Our definition also separates the usual concept of 
persistence into the two questions of whether trajectories are bounded and whether w-limit 
points exist on the boundary. Both properties are conjectured to hold for weakly reversible 
and, more generally, endotactic (see Definition I4.9P reaction networks. When trajectories 
are bounded, our definition of persistence is equivalent to the standard one. 

Remark 2.10. Suppose a manifold M has a metric d, and V is a compact subset of M. 
A differential inclusion X is repelled by V if for every d± € M>o, there exists c?2 € M>o such 
that d(f(I), V) > c?2 for every trajectory / : / — > M of X satisfying d(/(inf I), V) >d\. 

Definition 2.11. Differential inclusions that are persistent relative to the boundary dM are 
simply called persistent, and similarly for permanent and repelled. A collection of differential 
inclusions, possibly on a family of different manifolds with corners, is persistent, permanent, 
or repelled if each differential inclusion in the collection has the respective property. 
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Remark 2.12. Our use of the term "repelled" is related to other usages. Anderson and Shiu 
defined a boundary face to have a "repelling neighborhood" if the description in Remark l2.10l 
holds: if there is a neighborhood of the face such that whenever a trajectory is in that 
neighborhood, it can get no closer to the face while remaining in that neighborhood [2]. 
A different definition is that of a "repelling face" which Banaji and Mierczynski used to 
describe a boundary face for which any trajectory that begins in that face immediately exits 
the face into the interior of the relevant invariant set [6]. 

Remark 2.13. Suppose a differential inclusion X is a, subset of a persistent differential 
inclusion Y. Then X must be persistent, since its trajectories are all automatically trajec- 
tories of Y. More generally, consider properties P of differential inclusions X of the form 
u P(f) holds for all trajectories / of X." If property P is true for a differential inclusion Y, 
and a differential inclusion Z factors through Y, then property P is true for Z as well. 

3 Vertexical families of differential inclusions 

A mass-action kinetics system is naturally defined on the positive orthant M> correspond- 
ing to the space of concentrations of species. In this section we work instead with open 
hypercubes (0, l) s , and not directly with positive orthants themselves. 

To justify this choice, first we argue that little is lost by working with hypercubes. Open 
hypercubes (0, 1) are diffeomorphic to positive orthants R> , so differential inclusions can 
be transferred from one space to the other by fixing a diffeomorphism and using its Jacobian 
(see Section HT2]) . Additionally, properties such as persistence are defined topologically on 
the tangent bundle — and therefore invariant under diffeomorphism — so they can be analyzed 
on either space. 

Another advantage of working with open hypercubes is that they have natural "cubical" 
compactifications [0, l] s that appear to be optimally relevant in the context of mass-action 
kinetics. Alternatively, we could have achieved a cubical compactification by considering 
the hypercubes [0, oo] s . Nevertheless, there is a stylistic advantage to working with the 
hypercubes [0, l] s : we can treat symmetrically the cases where a species concentration 
goes to infinity or to zero. This makes some of our definitions more transparent, and the 
structural induction becomes cleaner. 

Definition 3.1. For any finite nonempty set S, let T>$ be the set of all differential inclusions 
on the open hypercube (0, l) s . Fix a collection S of finite nonempty sets. If X$ is a set of 
differential inclusions on (0, l) s for each S € S, then the collection X = {X$ C T>s}ses of 
sets X$ is a family of differential inclusions on open hypercubes indexed by S. 

3.1 Definitions concerning hypercubes 

The definition of vertexical families requires some preliminary notation on hypercubes. 
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Notation 3.2. Let S be a finite nonempty set. 

1. For every i G S, let G IR 5, be the map S — > R that sends i to 1 and S 1 \ {i} to 0. 

2. For each subset £7 C S and vertex x G {0, 1} of the hypercube [0, l] s , let 

F[/(x) = (x + span{ei | z 6 U}) D [0, 1] S 
denote the face of the hypercube [0, l] s along U at vertex x. 



3. For p = J2i£sPi e 'i ^ \p\ = y^iesPl denote its Euclidean norm. 

4. For subsets P,Q C M 5 , denote by d(P,<2) = inf { \p - q\ | p G P and g £ Q} the 
distance between them. 

In the context of reaction networks, S indexes the set of reacting chemical species; a 
chemical complex, being a linear combination of these species, is therefore viewed as a vector 
in R s (see Definition 14. ip . which has preferred basis vectors for i G S. 

Definition 3.3. Fix a finite nonempty set S. For a face F of [0, l] s and a real number 
rj > 1/2, the centered shrinking 

r]F = {x G F | d(x,<9F) > (1 - r{)/2} 

of F is the set of points in F whose distance from the boundary dF is at least (1 — r/)/2. 

Example 3.4. A centered shrinking of the rightmost face of the 3-cube looks as follows, 











; nF ; 








F/S 



where the inner shaded square is centered in F and has side length r/ times that of F . 

Definition 3.5. Fix a finite nonempty set S and a subset U C S. Let P C [0, l] 5 , and 
suppose that e G (0, 1/2) C M. The e-pile of the subset P along U is the set 

pile(P, e, £7) := jx + ^ x G P and - e < £i < e for all i G Uj fl [0, l] 5 . 
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Definition 3.6. Fix a finite nonempty set S and a proper face F of the hypercube [0, 1] 
containing a vertex x. Let U C S be such that i 7 = Fs\u( x )- For a real number e € (0, 1/2), 
the e-block F £ is pile ((1 — 2e)F, e, U) , the e-pile along U of the centered shrinking (1 — 2e)F. 

Example 3.7. If e = (1 — f?)/2 in Example 13.41 the block F e looks like the following: 



Note that the face F is orthogonal to the basis vector indexed by U, and the thickness e of 
the block F e equals its distance from the edges of F. The vertex x in Definition 13.61 could 
equally well be any of the four vertices of F. 

Remark 3.8. The block F £ is a closed subset of [0, l] s . Such sets are closely related to sets 
that Pantea p2] denoted by K £ , which he used for the purpose of projecting trajectories, 
as we too do in the current paper. 

Notation 3.9. Let / : U — > S be a map of finite sets, and view R s as functions S — > R. 
Denote by ttj : — > M, u the linear projection that sends v £ K s to v o f 6 MP . When 
U Q S and / is inclusion, we write ttjj for the projection map instead of 7ry. 

Remark 3.10. If U C 5 is any subset, then 

1. projection is surjective on the open hypercube: ^((0, l) s ) = (0,1)^, and 

2. 7T£/((1 — 2e)Fs\u{x)) is a vertex of the hypercube [0, l] u . 

Example 3.11. The projection ttjj in Example 13. 71 collapses the cube to the horizontal edge 



that is [0, 1] = [0, 1] . The projection takes F as well as the subset rjF C F to the indicated 
right-hand vertex of the interval. 




U 
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3.2 Definition of vertexical family and main result 

The heuristic description of a vertexical family of differential inclusions begins by considering 
a trajectory of a differential inclusion in the family. Suppose the trajectory remains near a 
face of the hypercube. The vertexical condition requires that the image of the trajectory 
under the projection map collapsing that face be the trajectory of a lower-dimensional 
differential inclusion in the family, up to time reparametrization. 

Definition 3.12 (Vertexical family). Let <S be the set of all finite nonempty subsets of the 
positive integers Z>i. A family X = {Xs}ses of differential inclusions on open hypercubes 
indexed by S is vertexical if for each 

• set S G S, 

• differential inclusion X C T(0, l) s in X$, 

• proper nonempty subset U C S, and 

• face F = F s \ v (x) of [0, l] s , 

there exists s > such that for every trajectory / : / — > [0, l] s of X with image in the 
block F e , there exist 

• a differential inclusion Y G Xjj, 

• a trajectory g : J — s> [0, l] u of Y, and 

• an order-preserving continuous map a : I —> J 

such that 7Tj/ o / = g o a. 

Examples of vertexical families of differential inclusions include those arising from re- 
versible, weakly reversible, endotactic, or strongly endotactic chemical reaction networks 
(Definitions 14.61 and 15. 15|) ; this is the content of Corollary 15.201 the goal of Sections H] and [5j 

Here is our main theorem on abstract vertexical families. 

Theorem 3.13. Fix a vertexical family X = {Xs}seS ° n open hypercubes indexed by the 
set S of all finite nonempty subsets of the positive integers Z>i. If for every set S € S, 
every differential inclusion X G X$ is repelled by the vertex set {0, 1} S of its hypercube, 
then every such X is repelled by, and hence persistent with respect to, the boundary <9[0, 1] . 

Proof. Fix S G S. If S has cardinality 1, then <9[0, l] s repels X$ by hypothesis. This 
constitutes the base case for an induction on the cardinality of S that allows us to assume 
Xjj is repelled by the boundary d[0, l] u for all proper nonempty sets U C S, with the goal 
of proving the result for a given differential inclusion X G X$. 

What we actually show is that for every proper, positive-dimensional face F = Fg\jj(x) 
of [0, l] s , if X is repelled by the boundary dF of the face then X is repelled by F itself. 
Given the hypothesis that X is repelled by the vertex set {0, 1} S , we deduce by induction 
on the dimension of F that every proper face of [0, 1] is repelling. Since the hypercube 
[0, l] s has only finitely many faces, its entire boundary is therefore repelling, as desired. 
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Suppose U C S is a proper nonempty subset, and let F = F s \u(x) be a corresponding 
face of the hypercube [0, 1] , parallel to the orthogonal complement of U. Assume that X 
is repelled by dF. Let /:/—>• (0, l) 5 be a trajectory of X, and let d\ = d(/(inf I),F) 
be the distance to F from the initial point of the trajectory. By Remark 12.101 it is enough 
to exhibit e > 0, depending only on d\, so that the trajectory remains at distance greater 
than e from F; that is, d(/(J),F) > e. 

Since X is repelled by dF, Remark 12.101 provides d 2 > such that d(f(I),dF) > d 2 . 
Decreasing c?2 if necessary, assume that d 2 < d\. Consider the block F d2 j 2 of F. If the 
image of the trajectory /(J) does not intersect the block F d2 / 2 , then by defining e = d 2 /2 it 
follows that d(/(I), F) > e, and we are done. Therefore, we can and do assume that I' C I 
is a nonempty maximal subinterval such that /(/') C F d2 / 2 , and let t = /(inf 7') denote the 
corresponding initial point, as in the following illustration. 



d 2 /2 

H — H 



5\J7 




Note that d(t, i 7 ) = d 2 /2. Indeed, by maximality of the interval /' the point t lies on the 
boundary of F d2 / 2 , and the only boundary face of F d2 / 2 that intersects the interior (0, l) s 
of the hypercube without also being contained in the ^-neighborhood of the boundary dF 
has constant distance d 2 /2 from F; see the following illustration. 



S\U 




u 
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Because X is a vertexical family, there exist a differential inclusion Y in Xjj, a trajectory 
g : J — > [0,1] U of y, and an order-preserving continuous map a : I' J such that 
g o a = ttjj o / with domain /'. It now suffices to show that there exists e > 0, depending 
only on d^ (which depended only on d\\ see Remark |2. lQj) such that d(g( J), ttu{F)) > e, 
for this claim implies that d (/(/), i 7 ) > e, as desired. 

To prove the claim, note by hypothesis that Y is repelled by the vertex iru(F) of [0, 1]^. 
Hence there exists e > such that all trajectories starting at distance g?2/2 away from the 
vertex ttjj(F) — and in particular, the trajectory g, because i has distance precisely g?2/2 
from F — never get closer than e to the vertex ttu(F). This choice of e depends only on the 
distance ^2/2, not on any aspect of the particular trajectory /. □ 

Remark 3.14. In the statement of Theorem 13.131 if the goal is to prove that a particular 
differential inclusion X £ X$ is repelled by (or persistent with respect to) <9[0, l] s , then it 
is enough to assume a slightly weaker hypothesis, namely that (i) X itself is repelled by 
the vertices {0, l} 5 of the hypercube, and (ii) the lower-dimensional set Xjj for each proper 
nonempty subset U C S is repelled by the vertices {0, 1}^ of its corresponding hypercube. 

4 Reaction network theory 

In this section, we define reaction networks and their properties. Our networks are more 
general than usual for chemical reaction network theory |1CH [T5] . 

Definition 4.1. Write J = {(a, b) \ < a < b < 00} for the set of all open subintervals 
of R>o and 3bdd = {[°> b] \ < a < b < 00} for the set of compact subintervals of M>o- 

1. A reaction network (S, C, TV) is a triple of finite sets: a set S of species, a set C C IR 5 of 
complexes, and a set 1Z C C x C of reactions. The reaction graph is the directed graph 
(C, TV) whose vertices are the complexes and whose directed edges are the reactions. 

2. A tempered reaction network is a reaction network (S, C, 1Z) together with a tempering 
K : 1Z — ^ 3bdd °f its reactions by nonempty compact positive intervals. 

3. A meted reaction network is a reaction network (S,C, TZ) together with an allotment 
fi : S — > 3 of an open positive interval to each species s £ S. The domain hypercube 
of a meted reaction network is the open hypercube D^ 1 = rises' M s ) != ^>o- 

A reaction r = (y,y') £ TZ, also written y — > y', has reactant y = reactant(r) G W s , product 
y' = product(r) £ M s , and reaction vector 

flux(r) = product(r) — reactant (r). 

The reaction diagram is the realization (C, 7£) — > M of the reaction graph that takes each 
reaction r £ 7£ to flux(r). A linkage class is a connected component of the reaction graph. 
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Remark 4.2. To explain the motivation for tempered networks, recall that a chemical re- 
action network gives rise to a dynamical system by way of reaction rates. For biochemical 
reaction networks, one is typically ignorant of the precise values for the rates. This oc- 
curs both because of incomplete information, and because of molecular and environmental 
variability. One way to model this ignorance is to allow reaction rates k(t) to be time- 
dependent, as long as they are uniformly bounded away from and oo. Craciun, Nazarov, 
and Pantea called such systems k- variable [7]. In a similar spirit, we have chosen to work 
with differential inclusions, allowing n(r) to take on every value from an appropriate interval. 

Remark 4.3. Intuitively, one thinks of a meted reaction network as a promise that concen- 
trations of species remain within the allotted intervals. The usual chemical reaction network 
theory setting would define /x as allotting the entire positive ray R>o to each species, so 
that the domain hypercube is the entire positive orthant D M = M> . 

Remark 4.4. The mathematical motivation prompting tempered and meted reaction net- 
works is to ensure that projections of trajectories "stay in the family" . Projections forget the 
exact concentrations of eliminated species. Absorbing the effect of pre-projection dynamics 
into post-projection dynamics requires a guarantee that the projected species concentrations 
never leave a certain sufficiently large interval. Therefore, in post-projection dynamics, the 
reaction rates remain within appropriately enlarged intervals. Tempered (interval-valued) 
reaction rates and meted reaction networks provide the extra flexibility for this construction. 

This intuition is made precise in Section [5j interval-valued rates allow us to define 
an appropriate domain category on which mass-action kinetics becomes functorial (Theo- 
rem ET7|). As a consequence, projective classes of reaction network (Definition 15. ip give 
rise to families of differential inclusions that are vertexical. 

Remark 4.5. The chemical reaction network theory literature usually imposes the following 
requirements for a reaction network: 

• each complex takes part in some reaction: for all y € C there exists y' € C such that 
(y, y') G 1Z or (y 1 , y) £ K, and 

• no reaction is trivial: (y, y) £ 1Z for all y € C. 

Definition 14.11 does not impose these conditions; in other words, our reaction graphs may 
include isolated vertices or self-loops. We drop these conditions to ensure that the projection 
of a network — obtained by removing certain species — remains a network under our definition 
even if some reactions become trivial (see Definition 15. 11 1). In addition, we allow arbitrary 
real complexes y € K s , so our setting is more general than that of usual chemical reaction 
networks, whose complexes y € Z> are nonnegative integer combinations of species, as in 
the following definition. 

Definition 4.6. A reaction network (S,C,1Z) is 
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1. integer if C C 7L S . 

2. chemical if C C Z> . 

3. reversible if the reaction graph of the network is undirected: a reaction {y, y') lies in 7£ 
if and only if its reverse reaction (y', y) also lies in 1Z. 

4. strongly connected if the reaction graph of the network is strongly connected; that is, 
if the reaction graph contains a directed path between each pair of complexes. 

5. weakly reversible if every linkage class of the network is strongly connected. 

Note that a network is strongly connected if and only if it is weakly reversible and has 
only one linkage class. 

Definition 4.7. The stoichiometric subspace H of a network is the span of its reaction 
vectors. For a positive vector xq € R>oj the invariant polyhedron of xq is the polyhedron 

V = (x + H ) n R§ . 

This polyhedron is also referred to as the stoichiometric compatibility class in the chemical 
reaction network theory literature. 

Definition 4.8. The standard basis of R 5 indexed by S defines a canonical inner product 
(•, •) on R 5 with respect to which the standard basis is orthonormal. Let w € R s . 

1. The vector w defines a preorder on R , denoted by < w , in which 

V<wV ^ (w,y) < (w,y'). 
Write y < w y' if (w,y) < {w,y'}. 

2. For a finite subset Y C R s , denote by init„, (Y) the set of ^^-minimal elements of Y: 

imt w (Y) = {y € Y \ (w,y) < (w,y') for all y' € Y}. 

3. For a chemical reaction network (S,C,1Z), the set 1Z W C TZ of w-essential reactions 
consists of those whose reaction vectors are not orthogonal to w. 

K w = {r eH\ (w,flux(r)) + 0}. 

4. The w-support supp^ (S, C, 1Z) of the network is the set of vectors that are <^-minimal 
among reactants of u;-essential reactions: 

supp w (S,C,TZ) = mit w (reactant(7^u,)) . 
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Definition 4.9. Fix a reaction network (S,C,H). 

1. The network (S,C,7Z) is w- endotactic for some w € R if 

(ti),flux(r)) > 

for all unessential reactions r G 7^ such that reactant(r) € supp^, (S,C,1Z). 

2. The network (S, C, TV) is W -endotactic for a subset W C R 5 if (S, C, TV) is ui-endotactic 
for all vectors w (zW. 

3. The network (S,C,TV) is endotactic if it is R 5 -endotactic. 

4. (S,C,TV) is strongly endotactic if it is endotactic and for every vector u> that is not 
orthogonal to the stoichiometric subspace of (S,C,1Z), there exists a reaction y — > y' 
in 7£ such that 

(i) y < w y', and 

(ii) y is ^-minimal among all reactants in (S,C,TV): y S init^ (reactant (TV) ) . 

Remark 4.10. Endotactic chemical reaction networks, which generalize weakly reversible 
networks, were introduced by Craciun, Nazarov, and Pantea [3 §4]. Our definition is 
slightly more general still, because we do not require the reaction networks to be chemical 
(Definition 14. 9p . Strongly endotactic reaction networks are new; they give rise to strong 
results concerning permanence using our techniques; see Theorem 16.51 

Remark 4.11. For the geometric intuition behind Definition 14. 9| imagine a hyperplane 
normal to w that is sweeping across the reaction diagram in R from infinity in direction w. 
As this hyperplane sweeps, it stops when it reaches the reactant y of a reaction y — > y' that 
is not perpendicular to w. If all such reactions do not point into the halfspace already swept 
by the hyperplane — that is, all such reactions have product y' outside of the open swept 
halfspace — then the network is u>-endotactic. Equivalently, the network is endotactic if no 
such reaction makes an obtuse angle with w. Illustrations can be found in [7]. 

As for strongly endotactic networks, the sweeping hyperplane now stops when it touches 
the reactant of any reaction, whether or not it is perpendicular to w. Again we require that 
the products of all such reactions lie outside of the open swept halfspace, and in addition 
at least one of these reactions is not perpendicular to w. If this condition is satisfied for 
all vectors w not orthogonal to the stoichiometric subspace, then the network is strongly 
endotactic. 

Example 4.12. Here we follow the usual convention of depicting a network by its reaction 
graph or reaction diagram and writing a complex as, for example, 2 A + B rather than 
y = (2, 1). The Lotka-Volterra reaction network, consisting of the three reactions 

A -> 2A A + B -> 2B B -> 0, 
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is not endotactic. Reversing all three reactions yields the network 

2A->A 2B->A + B (1) 
which is strongly endotactic, as can be verified from its reaction diagram: 



B 




A 



Example 4.13. Every weakly reversible reaction network is endotactic [7]. However, even 
a reversible reaction network may fail to be strongly endotactic, as in the following example 
of a pair of reversible reactions. For w = (0, 1), the < M -minimal reactant complexes are at 
the bottom, but neither of the corresponding reactions is perpendicular to w. 



B 




A 



5 Functoriality of mass-action kinetics 

In this section, we assign a differential inclusion to each meted tempered network (Defini- 
tion [5JJ)]), generalizing the usual mass-action kinetics ODE system in the chemical reaction 
network theory literature |1CH 1 15 j . Viewed appropriately, this assignment of differential 
inclusions gives mass-action kinetics a functorial nature (Theorem 15. lTj) . Consequently, 
projective classes of reaction networks (Definition 15. ip give rise to vertexical families of 
differential inclusions (Theorem 15. 19j) . In particular, the reversible, weakly reversible, en- 
dotactic, and strongly endotactic chemical reaction networks defined in the previous section 
all give rise to vertexical families (Corollary 15.201) . 

5.1 Categorical definitions 

Definition 5.1. Recall, from Definition 13. 9( the projection try : R s — > M. u for U C S, and 

denote by tt^ 2 = -k v x ttu : R s x R s ->■ M u x M u the product of ttu with itself. 

1. For a reaction network (S, C, 1Z), and a nonempty subset U C S of species, the reduced 
reaction network is the reaction network ttu(S,C,1Z) = (U, ttjj(C), 7r^ 2 (7£)) . 
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2. A property P of reaction networks is projective if for all finite nonempty sets S, 
reaction networks (S,C,TZ), and nonempty subsets U C S, if (S,C,1Z) has property P 
then the reduced reaction network iru(S,C,7Z) has property P. 

3. The set of all reaction networks with a given projective property is a projective class. 

Remark 5.2. The reduced reaction network iru(S, C, 1Z) is obtained from the reaction 
network (S, C, 1Z) by deleting all species outside of U. This concept was defined by Anderson 
[H §3.2], who required that any trivial reactions be removed from the reduced reaction set 
(ttu x ttjj)(TZ). In contrast, we allow trivial reactions (Remark I4.5|) . Another related notion 
in the context of reversible reactions is that of "reduced event-system" introduced in pQ. 

Example 5.3. For U = {A}, the reduced network 

2A -> A <- <- 

of the network ([1]) in Example 14. 121 is obtained by removing species B. The reduced network 
is strongly endotactic, as is the original network ([I]). 

The implication in Example 15.31 holds in general for strongly endotactic networks. 

Lemma 5.4. The classes of integer, chemical, reversible, strongly connected, weakly re- 
versible, endotactic, or strongly endotactic reaction networks are projective. Further, if P\ 
and Pi are projective properties, then so are the conjunction P\t\Pi and disjunction Pi V Pi- 
Proof. Projectivity holds for integer and chemical networks because projection preserves 
integrality and nonnegativity of points in W s . Projectivity holds for reversible, strongly 
connected, and weakly reversible networks because these conditions depend only on the 
reaction graph, on whose vertices and edges projection is surjective. 

Next, consider an endotactic network with species set S and the reduced network arising 
from a nonempty subset U C S. Take any vector w € MF . For any reduced reaction nu(r), 
where r is a reaction in the original network, 

(it;,flux(7rtf(r))) = (K0),flux(r)>. (2) 

Thus, the unessential reactions of the reduced network are the projections under 7r^ 2 of the 
(w, 0)-essential reactions of the original network, where we write (w,0) & M. u x R s \ u . Simi- 
larly, the u>-support of the reduced network is the projection under ttjj of the (w, 0)-support 
of the original network. So, if nu(r) is a unessential reaction of the reduced network with 
reactant in the u-support, then the original reaction r is a (u, 0)-essential reaction of the 
original network with reactant in the (w, 0)-support. By ([2]) and the definition of endotactic, 
(u>, flux(7T[/(r))) = ((u, 0), flux(r)) > 0. Hence the reduced network is endotactic. 

Next, let H denote the stoichiometric subspace of a strongly endotactic network, so 
ttjj{H) is the stoichiometric subspace of the reduced network. Take a vector w € M. u that is 
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not orthogonal to ttu(H). We need only show that there exists a reaction iru(y) — > iru (y 1 ), 
where y — > y' is a reaction in the original network, such that iru(y) < w ^u{y') an d t^u{v) 
is ^^-minimal among all reactant vectors in the reduced network. Again consider (w, 0) G 
M. u x M 5 ^. As w is not orthogonal to iru(H), it follows that (w,0) is not orthogonal to H, 
and the preorder < w on the reduced reactant complexes is the projection under -kjj of the 
preorder <( w ,o) on the original reactant complexes. Since the original network is strongly 
endotactic, there is a reaction y — > y' in the original network with y <( w ,o) v' such that y is 
<( U)! o)- m uh ma l among all reactant vectors. This reaction achieves our requirements. 

The claim about conjunctions and disjunctions follows formally by Definition l5.1l2i □ 

Notation 5.5. Let I,JQ M>o be two intervals, and S a finite nonempty set. 

1. Define / • J = {i ■ j \ i G /, j G J} C K>o and Q s£S I s pointwise. 

2. For n G Z>i define I n recursively as / • I n ~ 1 with I 1 = I. 

3. I x J and IlseS-^ denote Cartesian products of intervals, as usual. 

Definition 5.6. Let S be a finite nonempty set, and consider a function fi : S — > 3 to the 
set of open positive intervals. A nonempty subset U C S is fi-projectable if < inf n(s) and 
sup /i(s) < oo for all s G S\ U; that is, the left and right endpoints of the intervals /x(s) are 
bounded away from and oo for those s outside of U. 

Remark 5.7. The condition of Definition 15.61 is on the complement S \ U because those 
are the species removed in projecting to U, and so it is those species that must be bounded 
away from and oo. The set S itself is trivially ^i-projectable, for all \l : S — >■ 3. 

Intuitively, meted tempered networks form a category whose morphisms are projections, 
where the projection pu from one object N to another another corresponds to substituting 
intervals from the allotment of N in place of a set S \ U of forgotten species. 

Definition 5.8. The category M of meted tempered reaction networks with projections is 
given by the following data. 

1. Objects: each is an intervaled tempered reaction network N, specified by a reaction 
network (S, C, 1Z) along with a tempering k : TZ — > 3bdd an d an allotment /i : S — > 3. 

2. Morphisms: pu : N ->• N' if 

• the network of N' is (S',C',1Z') = irjj(S,C,lZ) for a ^-projectable subset U C S; 

• the tempering of N' is k' : 7T[/(r) i-> n(r) ■ ^ s ^reactant(r) s ^ w j iere exponent 

ses\u 

on fj,(s) is the component indexed by s in the vector reactant(r); and 

• the allotment of N' is \s! = fi\jj, obtained by restricting the allotment of N to U. 
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Remark 5.9. Composition in M is well-defined because projecting to U C S and then 
projecting to V C [7 is the same as projecting directly to V. 

In our conception, mass-action kinetics produces from each meted tempered reaction 
network a differential inclusion on its domain hypercube. Our key observation along these 
lines, Theorem 15.171 states that mass-action kinetics is a functor, with domain category AT 
and codomain category as follows. 

Definition 5.10. The category T>Z of differential inclusions is given by the following data. 

1. Objects: each is a choice of manifold with corners and a differential inclusion on it. 

2. Morphisms: a morphism from X C TM to Y C TN is a continuous map k : M — >• N 
such that for each trajectory / : I — )■ M of X, there is a trajectory <? : J — >■ iV of Y 
and an order-preserving continuous map a : I — > J satisfying 

ko f = goa. (3) 

Lemma 5.11. Composition of continuous maps induces a well-defined composition on DZ. 
Specifically, assume Xj C TMj /or j = 1,2,3 are differential inclusions, with morphisms 
k\2 : M\ — > M2 and k^z : M2 — > -M3 in DZ. If fi : I\ — > M\ is a trajectory of X\, then there 
is a trajectory fa : I3 — > M3 0/ X3 and an order-preserving continuous map 013 : I\ — > 1-$ 
such that the composite continuous map k\^ = k<iz o kyi satisfies k\-& f\ = fz° 013- 

Proof. Given fi, since k±2 is a morphism in DZ, there is a trajectory fi : I2 — > M2 of X2 
and a continuous order-preserving map a±2 : I± —¥ I2 such that k\2 /1 = /2 «12- For the 
desired trajectory fa : J3 — > M3 of X3 use the one afforded by virtue of k23 being a morphism 
in T>Z, given fa, which comes with a continuous order-preserving map 023 : fa — > H such 
that fc 2 3 fa = fa «23- Set ai3 = a 2 3 «12- Then 

kl3 fl = k 2 3 ^12 fl 
= k 2 3 o fa° Ct\2 
= fa° «23 "12 
= fa "13, 

as desired. □ 

Remark 5.12. Compare the notion of morphism in T>Z (Definition I5.10p with that of 
vertexical family (Definition 13. 12p . Equation ([3]) also occurs in Definition 13. 12[ with k being 
a particular type of continuous map ttu- However, Definition 15.101 asks for a global map k, 
whereas the maps in Definition 13.121 are required only locally, on blocks of faces. 
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Remark 5.13. One motivation for defining the category of differential inclusions this way 
comes from the dynamical systems concept of topological equivalence |16j . which identifies 
two phase portraits as qualitatively the same, even if the details of the dynamics may differ. 
The isomorphisms in our category T>X correspond exactly to topological equivalence. 

For this reason, morphisms between differential inclusions may also be called topological 
morphisms. Intuitively, if a topological morphism is a monomorphism, then its target differ- 
ential inclusion qualitatively simulates the domain differential inclusion. Maps that are not 
monomorphisms can of course result in the loss of information, in general. The categorical 
message of Theorem 13.131 is that it is sometimes possible to piece together many "lossy" 
maps on the same domain to regain substantial information about the domain dynamics. 

Another concept from dynamical systems, topological congruence [16], is a stronger 
notion than topological equivalence that disallows time reparameterization. This motivates 
looking at a subcategory T>X\ of T>X in which the order-preserving map a in Definition 15 .101 
is required to be the identity map; the definition follows. 

Definition 5.14. The category T>X\ of differential inclusions with topological congruence 
morphisms is the subcategory of T>X with the following data. 

1. Objects: the same objects as in T>X. 

2. Morphisms: a morphism from X C TM to Y C TN is a continuous map k : M — > N 
such that k o f is a trajectory of Y whenever / : I — s> M is a trajectory of X. 

The proof of Lemma 15.111 makes it plain that any composition in T>X of morphisms 
in T>X\ is a morphism in because every reparameterization map a™ in that proof can 
be taken to be the identity map on I\. 

Definition 5.15. The mass-action differential inclusion of a meted tempered reaction 
network N, specified by a reaction network (S,C,1Z) with tempering k and allotment fi, is 
the differential inclusion on the domain hypercube D M C whose fiber over x € D M is 

{ J2 k r flux(r):c rcactant ( r ) | kr G «(r) for all r € ft} C R s = T X W. 

The mass- action functor 9Jt : M — > T>X\ is defined on 

1. objects N by setting 9JT(iV) equal to the mass-action differential inclusion of N, and on 

2. morphisms pjj : N — > N' by setting d3l(pu) equal to the projection irjj '■ — > D^'- 

Remark 5.16. The usual mass-action kinetics ODEs are the special cases where the tem- 
pering k assigns not an interval of positive length but a point — the reaction rate constant — to 
each reaction, and the domain hypercube is the entire positive orthant: = M^ . 
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5.2 Functor ial results and consequences 



Theorem 5.17. The mass-action functor 9Jt is a functor from the category J\f of meted 
tempered reaction networks with projection morphisms to the category T>X\ of differential 
inclusions with topological congruence morphisms. 

Proof. The content of the statement is twofold: first, that iJJl(pu) = i^u m Definition 15. 151 2 
indeed defines a morphism ^Xft(N) — > yyt(N'), and second, that 50? preserves identity mor- 
phisms as well as compositions. The second is straightforward (projections are sent to 
projections). For the first, it suffices to observe that for any nonempty U C S, the projec- 
tion ttjj o f of a trajectory / of 9Jt(N) is a trajectory of 9Jt(pu(N)) by definition of pjj. □ 

Vertexical families are defined in terms of differential inclusions on unit hypercubes 
(0, l) s . In contrast, the mass-action functor produces differential inclusions on positive 
orthants M^ . To translate back and forth, our next definition fixes a smooth, order- 
preserving diffeomorphism R>o — > (0,1). The actual choice is irrelevant for our purposes, 
but if it helps, the reader may consider the function x t— >■ x/(l + x). 

Definition 5.18. Fix a smooth, order-preserving diffeomorphism I : R>o — s- (0, 1). For 
every nonempty finite set S, let £ s : R$ -> (0, l) 5 , with derivative d£ s : TR§, T(0, l) 5 . 

1. The pushforward of a differential inclusion X on R> under I is the differential inclu- 
sion dl s {X) on (0,l) s . 

2. For a meted tempered reaction network N defined on a species set S with allotment /i, 
the differential inclusion dJl i (N) is the pushforward of the mass-action differential 
inclusion VJt(N), considered as a differential inclusion on the image ^(□ M ) C (0,1) 5 
of the domain hypercube under I. 

The following is the main result of this section. 

Theorem 5.19. Let P be a projective property of reaction networks. If Tp is the class 
of all meted tempered reaction networks whose underlying reaction networks have property 
P and whose domain hypercubes equal the entire (open) positive orthant, then Wl^^p) = 
{9Jl(N) | N € J~p} is a vertexical family of differential inclusions on open hypercubes. 

Proof. By definition, SJT^J-p) is a family of differential inclusions on open hypercubes. 
Hence we need only prove 9Jt e (J r p) is vertexical. Consider a meted tempered reaction net- 
work N € J-p, specified by a reaction network (S, C,TZ) with tempering k and allotment \i. 
Fix a proper nonempty subset U C S and a vertex x £ {0, 1} S of the hypercube (0, l) s . 
Denote by F = F s \u(x) the corresponding face of the hypercube. 
Fix < e < \ , and define // = jJ E v : S — >■ 3 by 





fi(s) if s € U 

£- l {e,l-e) if s eS\U 
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recalling that /x(s) = (0, oo) for all s € 5. Denote by N' the meted tempered reaction 
network N' that agrees with N except that the allotment of N' is // instead of /x. The 
underlying reaction network of N' still has property P, because ./V and N' have the same 
underlying reaction network. 

While N' itself no longer has the entire positive orthant M> as its domain hypercube, 
the set U is //-projectable (Definition 15.61) . and the reduced network pu(N') has domain 
hypercube K> - Since P is projective, pu(N') therefore lies in Tp. Consequently, by 
definition of vertexical family (Definition 13. 12H . it is enough to show that for any trajectory 
/ : / — s> F £ of 97t e (N) with image in the block F £ , the projection pjj o f = g is a trajectory 
of the mass-action differential inclusion 9Jl e (pu(N')). 

The trajectory / is also a trajectory of Wl e (N'), since F £ C I s ) by construction of \j! . 
The result is now deduced easily from functoriality of mass-action kinetics (Theorem 15. 17|) : 
the morphism ttjj in the category T)X\ from dJt(N') to 9Jt(pjj(N')) yields the pushforward 
morphism 9Jl e (N') — > yjl i (pu(N')) in T>T\, so the proof is complete by definition of mor- 
phisms in the category T>X\ (Definition 15.141 2). □ 

Recall that in propositional logic, a monotone (or monotonic) formula is one formed by 
the application of AND and OR operations only, without the use of NOT operations. 

Corollary 5.20. The classes of (monotone combinations of) integer, chemical, reversible, 
strongly connected, weakly reversible, endotactic, or strongly endotactic tempered reaction 
networks on positive orthants generate vertexical families of differential inclusions. 

Proof. Immediate from Theorem 15.191 and Lemma 15.41 □ 

Remark 5.21. The category N of networks with projections suffices for our purposes, 
but it would be more natural to allow tempered reaction networks with an associated 
domain set D C ]R^ that is not necessarily a Cartesian product of intervals. In addition, 
it is tempting to add to the category more morphisms, such as those corresponding to 
translation of the reaction diagram within R 5 , or scaling, rotation, arbitrary linear maps, 
graph homomorphisms, inversion z = 1/x, and so on. It is easy to verify that translation 
acts as a time-reparametrization. Hence, even allowing translations, mass-action kinetics 
remains a functor to T>X. This leads to the following problem, which we leave open. 

Question 5.22. What is the richest domain category for which mass-action kinetics re- 
mains a functor to the category T>T of differential inclusions? 

This question is important because a richer domain category would imply more ways 
of reducing the behavior of a network's mass-action kinetics to the behaviors of related 
networks. This could allow us to "program" (and analyze) instances of reaction dynamics 
in high dimensions as appropriate combinations of simpler reaction dynamics. 
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Remark 5.23. Consider a vertexical family of mass-action differential inclusions for which 
the one-dimensional differential inclusions in the family are known to be permanent. For 
instance, the differential inclusions arising from weakly reversible, endotactic, or strongly 
endotactic networks have this attribute. It is tempting to attempt to argue that such a 
family is permanent by the following induction: given a trajectory, all of its one-dimensional 
projections (which are also in the family, due to the functoriality of mass- action kinetics 
and projectivity of the relevant properties) are permanent, and hence the trajectory itself 
must be permanent. However, this does not work because of uniformity issues. It is true 
that for a given trajectory, there exists a compact set that it enters eventually. However, 
what we need is one compact set so that every trajectory eventually enters this set, and our 
inductive argument does not prove this. A successful argument about permanence requires 
additional structure; for example, see Remark 16.41 

6 Implications for persistence of mass-action systems 

One of the long-standing open problems of chemical reaction network theory is the global 
attractor conjecture concerning so-called "complex-balanced systems". Complex-balanced 
systems form a well-studied subclass of weakly reversible mass-action ODE systems that 
contain all so-called "detailed-balanced" systems and weakly reversible "deficiency zero" 
systems. Many properties about complex-balanced systems (as well as detailed-balanced 
systems and deficiency zero systems) were elucidated by Feinberg, Horn, and Jackson in 
the 1970s, and we provide only an overview here. (See any of the references [H [9] [TO], [15] 
for a definition of complex-balanced systems.) 

For complex-balanced systems, it is known that there is a unique steady state within the 
interior of each invariant polyhedron V. This steady state is called the Birch point in [8] due 
to the connection to Birch's Theorem in algebraic statistics. Moreover, a strict Lyapunov 
function exists for this point, so local asymptotic stability relative to V is guaranteed [101115] 
(see Remark 16. 4|) . An open question is whether all trajectories with an initial condition in 
the interior of V converge to the unique Birch point of V. The assertion that the answer is 
"yes" is the content of the following conjecture, which was stated first by Horn in 1974 [14] 
and given the name "Global Attractor Conjecture" by Craciun, et al. [8]. 

Conjecture 6.1 (Global Attractor Conjecture). For any complex-balanced mass-action 
system (S, C,1Z, k) and any strictly positive initial condition x° , the Birch point x € V := 
(x°+i?)n]Rf is a global attractor of the relative interior of the invariant polyhedron int('P). 

In 1989, Feinberg [11] conjectured that all mass-action ODE systems arising from weakly 
reversible reaction networks are persistent, and observed that the persistence of weakly 
reversible reaction networks would imply the global attractor conjecture. This reduced the 
problem of proving the global attractor conjecture to the problem of establishing persistence 
for weakly reversible mass-action systems. Our functoriality results reduce persistence to 
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bounding the dynamics away from vertices of the hypercube — so it suffices to ensure that 
all species remain bounded away from and oo in the positive orthant — at the price of 
considering tempered reaction networks. 

Corollary 6.2. Let J- be the class of all tempered reaction networks that are one or more of 
the following: integer, chemical, reversible, strongly connected, weakly reversible, endotactic, 
or strongly endotactic. If the mass-action differential inclusions on positive orthants of all 
reaction networks in T are repelled by vertices after pushing forward to open hypercubes by 
a smooth order-preserving diffeomorphism, then they are repelled by the entire boundary. 

Proof. Immediate from Theorem 13.131 and Corollary 15.201 □ 

Recently, Craciun, Nazarov, and Pantea [7] generalized Feinberg's persistence conjecture 
in the following three ways: the weakly reversible hypothesis is weakened to endotactic, 
fixed reaction rate constants are allowed to vary within bounded intervals (i.e., the network 
is tempered), and the conclusion of persistence is strengthened to permanence §4]. We 
present an extended version of their conjecture that drops the requirement that the reaction 
network be chemical. 

Conjecture 6.3 (Extended Permanence Conjecture). If J- is the family of meted tempered 
endotactic reaction networks whose domain hypercube is the entire positive orthant K> , 
then 9Jt e (J-) is permanent. 

Remark 6.4. Horn and Jackson [15] established that any complex-balanced mass-action 
ODE system N on a species set S admits the strict Lyapunov function 

9a(x) = y^Xi ( log — - 1 

where a G M> is a given complex-balanced steady state of N. Consequently, the differential 
inclusion 9Jl e (N) is vertex-repelling (see Definition 15. 18p . The proof, which is due in part to 
Anderson and Craciun et al. [3 [8], proceeds by showing that the Lyapunov function g a has 
a local maximum at the vertex as well as bounded level sets that do not allow trajectories 
to escape to infinity. 

If the class of complex-balanced systems were projective, then our arguments would 
have proved the global attractor conjecture. However, this is not the case: projections of 
complex-balanced systems need not be complex-balanced. Endotactic networks, however, 
do constitute a projective class, and, in our view, it is this property that allowed Craciun, 
Nazarov, and Pantea to make their projection-type arguments [7} 118] . Similarly, the prop- 
erty of having only one linkage class is projective: a network with only one linkage class 
maintains this property after reduction. This was used in Anderson's proof of the global 
attractor conjecture for networks possessing only one linkage class [3|. 
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Conjectures 16.11 and 16.31 both remain open; for an overview of recent progress on these 
problems, we refer the reader to the work of Anderson [3] §1.1]. However, using the results 
presented in the current paper, we prove the following in a subsequent paper |13j . 

Theorem 6.5. If J 7 is the family of meted tempered strongly endotactic reaction networks 
whose domain hypercube is the entire positive orthant M> > ^ en is permanent. 

The key step contributed by the results in the current paper is Corollary 16.21 in the case 
of strongly endotactic networks. The approach in [13] shows that for the mass-action differ- 
ential inclusions of strongly endotactic networks, outside a compact set the function g a {x) 
in Remark 16.41 continues to behave like a Lyapunov function. Theorem [63] is established by 
an argument along the lines suggested in Remark 16.41 
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